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ABSTRACT 

We present a multi- wavelength study of the Magehanic Stream (MS), a massive gaseous structure 
in the Local Group that is believed to represent material stripped from the Magellanic Clouds. We 
use uhraviolet, optical and radio data obtained with HST/COS, VLT/UVES, FUSE and CASS 
to study metal abundances and physical conditions in the Stream toward the quasar Fairall9. Line 
absorption in the MS is detected from a large number of metal ions and from molecular hydrogen (H2). 
The UVES data indicates a complex velocity-component structure with seven absorption components 
related to the MS. From the analysis of unsaturated Sii absorption, in combination with a detailed 
photoionization model, we obtain a surprisingly high a abundance in the Stream toward Fairall9 of 
[S/H] = —0.30 ± 0.04 (0.50 solar). This value is 5 times higher than what is found along other MS 
sightlines based on similar COS/UVES data sets. In contrast, the measured nitrogen abundance is 
found to be substantially lower ([N/H] = —1.15 ± 0.06), implying a very low [a/N] ratio of —0.85 
dex. The elements Si, Fe, Ni, Ti, Al are depleted into dust grains with depletion values in the range 
S « —0.5 to —1.6 dex. The H2 absorption in the Stream is accompanied by very narrow optical 
Nai/Caii absorption. Our analysis of these species suggests the presence of a cold (T ~ 100 K) 
and dense (njj ~ 5 cm~^) gas phase in the body of the Stream, possibly arranged in a widespread 
population of pc-scale clumps or filaments. From the analyis of C 11* absorption we infer a local 
electron density of rie < 0.05 cm~^ and a mean 0+ cooling rate of 7.3 x 10~^^ ergs"^. The substantial 
differences in the chemical composition of the Magellanic Stream toward Fairall 9 compared to other 
lines of sight, together with the presence of dust and molecular gas, point toward a complex enrichment 
history of the Magellanic Stream. We favour a scenario, in which the gas was locally enriched with a 
elements by massive stars in a star-forming region in the SMC or LMC and then was expelled and/or 
stripped from the stellar disk before the delayed nitrogen enrichment from intermediate-mass stars 
could set in. 
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Subject headings: ISM: abundances - Galaxy: halo - Galaxy: evolution - Magellanic Clouds - quasars: 
absorption lines 



1. INTRODUCTION 

The distribution of neutral and ionized gas in the cir- 
cumgalactic environment of galaxies is known to be an 
important indicator of the past and present evolution 
of galaxies. Both the infall of metal-poor gas from inter- 
galactic space and from satellite galaxies and the outflow 
of metal-rich gaseous material through galactic winds 
represent key phenomena that determine the spatial dis- 
tribution and the physical state of the circumgalactic gas 
around massive galaxies. 

From observations and theoretical studies, it is known 

^ Based on observations obtained with the NASA/ESA Hubble 
Space Telescope, which is operated by the Space Telescope Science 
Institute (STScI) for the Association of Universities for Research in 
Astronomy, Inc., under NASA contract NAS5D26555, and on ob- 
servations collected at the European Organisation for Astronomical 
Research in the Southern Hemisphere, Chile. 



that galaxy interactions between gas-rich galaxies can 
transport large amounts of neutral and ionized gas into 
the circumgalactic environment of galaxies. In the local 
Universe, the most massive of these extended tidal gas 
features can be observed in the 21 cm line of neutral hy- 
drogen (Hi). The most prominent nearby example of a 
tidal gas stream produced by the interaction of galaxies is 
the Magellanic Stream (MS), a massive (~ 10^-10^ Mg) 
stream of neutral and ionized gas in the outer halo of the 
Milky Way at a distance of ~ 50 — 60 kpc (e.g., Wan- 
nier & Wrixon 1972; Gardiner & Noguchi 1996; Weiner 
& Wilhams 1996; Putman et al. 2003; Briins et al. 2005; 
Fox et al. 2005, 2010; Koerwer 2009; Besla et al. 2007, 
2010, 2012). The Magehanic Stream spans over 200° 
on the sky (e.g., Nidever et al. 2010) and has a (mean) 
metallicity that is lower than that of the Milky Way, but 
comparable with the metallicity found in the SMC and 



p. Richter et al. 



LMC (0.1 - 0.5 solar; Lu et al. 1994; Gibson et al. 2000; 
Sembach et al. 2001; Fox et al. 2010, 2013). The MS also 
contains dust grains and diffuse molecular hydrogen (H2 ; 
Sembach et al. 2001; Richter et al. 2001a). A number 
of theoretical studies, including tidal models and ram- 
pressure stripping models, have been carried out to de- 
scribe the Stream's motion in the extended halo of the 
Milky Way and pinpoint its origin in one of the two Mag- 
ellanic Clouds (Gardiner & Noguchi 1996; Mastropictro 
et al. 2005; Connors ct al. 2006; Besla et al. 2010; Diaz & 
Bekki 2011). 

The origin and fate of the Magellanic Stream is closely 
related to the trajectories of LMC and SMC (e.g., Con- 
nors et al. 2004, 2006; Besla et al. 2007), and any realistic 
model of the MS thus needs to consider the dynamical 
and physical state of the Milky Way/Magellanic Clouds 
system as a whole (see also Bland- Hawthorn et al. 2007). 
While early tidal models have assumed that the Mag- 
ellanic Stream is a product from the tidal interaction 
between LMC and SMC as they periodically orbit the 
Milky Way (e.g., Gardiner & Noguchi 1996), more recent 
proper motion measurements of the Magellanic Clouds 
(MCs) (Kalhivayahl et al. 2006a, 2006b) indicate that 
the MCs may be on their first passage around the Milky 
Way. Some subsequent tidal models (Besla et al. 2010; 
Diaz & Bekki 2011) thus favour a first-infall scenario for 
the Magellanic Stream. Moreover, while many models 
(e.g., Connors et al. 2006) place the origin of the Stream's 
gaseous material in the SMC, other, more recent studies 
trace back at least part of the Stream's gaseous mate- 
rial in the LMC (e.g., Nidever et al.2008). The latter 
study also highlights the role of energetic blowouts from 
star-forming regions in the LMC for the formation of the 
Stream. Clearly, further theoretical studies and observa- 
tions are required to pinpoint the origin of the MS based 
on different (independent) methods. 

In the first paper in our series analyzing the chem- 
ical and physical conditions in the Magellanic Stream 
(Fox et al. 2013; hereafter Paper I), we have investi- 
gated MS absorption in the UV and optical along the 
hues of sight toward RBS 144, NGC 7714 PHL 2525, and 
HE 0056— 3622. In this paper we analyze the MS using 
UV and optical absorption-line spectra of the Seyfert 1 
galaxy Fairall9 (zom — 0.047). Located at / = 295.1 and 
b ~ —57.8 the Fairall9 sightline lies only 14.3° on the 
sky from the SMC. This sightline is the best-studied in 
absorption of all MS directions (Songaila 1981; York et 
al. 1982; Lu, Savage & Sembach 1994; Gibson et al. 2000; 
Richter et al. 2001a; Sembach et al. 2003), largely because 
the Fairall9 is bright in both the optical and the UV 
and the Stream's H i column in this direction is large 
(log N(Ri) « 20; see Gibson et al.2000). The high col- 
umn of neutral gas ensures that a wide range of low- 
ionization UV metal lines arc detectable in the Stream, 
and even molecular hydrogen was observed in the MS to- 
ward Fairall 9 data from the Far Ultraviolet Spectroscopic 
Explorer {FUSE; Richter et al. 2001a; Wakker 2006). Us- 
ing a spectrum of Fairall 9 obtained with the Goddard 
High Resolution Spectrograph (GHRS) onboard the Hub- 
ble Space Telescope (HST) together with Parkes 21 cm 
H I data Gibson et al. (2000) derived a metallicity of the 
Stream toward Fairall 9 of [S/H] = -0.55 ± 0.06i^:2l 
(~ 0.3 solar), which represented the most accurate metal- 



licity determination of the Stream from UV absorption 
line data at that time. This metallicity is consistent with 
either an SMC or LMC origin of the gas. A difficulty for 
constraining the origin the MS in one or the other Mag- 
ellanic Cloud arises from the fact that the gas in the 
Stream was stripped from its parent galaxy ~ 1 — 2 Gyr 
ago (e.g., Gardiner & Noguchi 1996; Connors et al. 2006; 
Nidever et al. 2008), but has not experienced any further 
metal enrichment since then, while the parent galaxy un- 
derwent further chemical evolution. The MS does not 
contain any massive stars (e.g., Mathewson et al. 1979), 
in contrast to the Magellanic Bridge (Irwin, Kunkel & 
Demers 1985). This aspect needs to be taken into ac- 
count when comparing metal abundances in the Stream 
with present-day LMC and SMC abundances. 

To increase the accuracy of the metallicity determina- 
tion of the Stream toward Fairall 9 and to obtain more 
detailed information on the chemical composition of the 
gas and dust in the Stream, more accurate spectral data 
are desirable. Because the Magellanic Stream is a mas- 
sive gas cloud with complex internal kinematics (e.g., 
Nidever et al.2008), data with high spectral resolution 
and a high signal-to-noise (S/N) ratio are required to 
fully resolve the Stream's velocity-component structure 
and to detect weak absorption features from the various 
metal ions that have their transitions in the ultraviolet 
(UV) and in the optical. As part of our ongoing project 
to study the properties of the Magellanic Stream in ab- 
sorption along multiple lines of sight (see also Paper I) we 
obtained high-resolution optical data of Fairall 9 from the 
Ultraviolet and Visible Echelle Spectrograph (UVES) in- 
stalled on the Very Large Telescope (VLT) and medium- 
resolution UV data from the Cosmic Origins Spectro- 
graph (COS) onboard the HST, both data sets providing 
absoption spectra with excellent S/N ratios. The combi- 
nation of these data sets, as desribed in this study, there- 
fore provides a particular promising strategy to study in 
great detail the chemical and physical conditions in the 
Magellanic Stream in this direction. 

This paper is organized as follows: in Sect. 2 we de- 
scribe the observations and the data reduction. The col- 
umn density measurements and the profile modeling are 
explained in Sect. 3. In Sect. 4 we derive chemical and 
physical properties of the gas in the MS. We discuss our 
results in Sect. 5. Finally, a summary our study is given 
in Sect. 6. 

2. OBSERVATIONS AND SPECTRAL ANALYSIS 

2.1. VLT/UVES observations 

FairaU9 was observed with the VLT/UVES spectro- 
graph (Dekker et al. 2010) in 2010 under ESO program 
ID085.C-0172(A)(PI: A. Fox). The observations were 
taken in Service Mode using Dichroic 1 in the 390-f 580 
setting, a 0.6" slit, and no rebinning. The observations 
were carried out under good seeing conditions (<0.8"). 
The raw data were reduced with the standard UVES 
pipeline, using calibration frames taken close in time to 
the corresponding science frames. The reduction steps 
involve subtraction of the bias level, inter-order back- 
ground, sky background, night sky emission lines, and 
cosmic ray hits. The frames were then flat-fielded, opti- 
mally extracted and merged. The wavelength scale was 
corrected for atmospheric dispersion and heliocentric ve- 
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locity and then placed into the local standard of rest 
(LSR) velocity frame. Multiple exposures on the same 
target were registered onto a common wavelength grid 
and then added. The final spectra have a very high 
spectral resolution of i? « 70 000 corresponding to a 
FWHM of 4.3 kms^^. They cover the wavelength range 
between 3300 and 6800 A. The S/N ratio per resolution 
element is 40 at 3500 A (Tin), 65 at 4000 A (Can), and 
83 at 6000 A (Nai). The UVES data thus provide much 
higher sensitivity and substantially higher spectral res- 
olution than previous optical measurements of Fairall 9 
(Songaila 1981). 

2.2. UST/COS observations 

Fairall 9 was observed with the HST/COS spectro- 
graph (Green et al. 2012) in 2012 under HST program 
ID 12604 (PI: A. Fox). A four-orbit visit provided a 
total of 5378 s of exposure time with the G130M/1291 
wavelength setting, and 6144 s with the G160M/1589 set- 
ting. With each grating, all four FP-POS positions were 
used to dither the position of the spectrum on the de- 
tector to reduce the fixed-pattern noise. The raw data 
were processed and combined with the CALCOS pipeline 
(v2.17.3). For the coaddition of the individual exposures 
we used interstellar absorption lines as wavelength refer- 
ence. The final, co-added spectra then were transformed 
into the LSR velocity frame. The COS data cover the 
UV wavelength range between 1131 and 1767 A. The 
spectra have a resolution of i? w 16 000 (FWHM « 19 
kms~^) and a pixel size of '^ 7 kms~^. The S/N ratio 
per resolution element is 37 at 1200 A and 26 at 1550 A. 

In order to minimize geocoronal emission, which con- 
taminates the absorption lines of Oi A1302, Sin A1304 
and Hi Lya in the velocity range —200 kms~^ ^^^lsr^ 
-f200 kms^^ during orbital daytime, we re-reduced the 
data with a night-only extraction. For this, data were ex- 
tracted from those time intervals when the Sun's altitude 
was less than 20 degrees. 

2.3. FUSE observations 

As part of our study we also re-analyze archival FUSE 
spectra of Fairall 9. These spectra were obtained in 2000 
under FUSE program ID PlOl (PI: K.R. Sembach) with 
a total exposure time of 34827 s. They show strong 
molecular hydrogen absorption arising in the MS, as pre- 
sented by Richter et al. (2001a). The FUSE spectra have 
a resolution of ~ 20 kms~^ (FWHM), and cover the 
wavelength range 912 — 1180 A. The raw data were re- 
duced using the CALFUSE pipeline v3.2.1. The individ- 
ual exposures were carefully coadded using interstellar 
lines as wavelength reference. Unfortunately, the S/N in 
the FUSE data is very low, which severely hampers the 
analysis of absorption lines (see also Richter et al. 2001a; 
Wakker 2006). From a detailed inspection of the indi- 
vidual spectra in the different segments we conclude that 
only the data from the LiF lA segment can be properly 
coadded without introducing large systematic errors in 
the flux distribution in the spectrum. In the coadded 
LiF lA data, that are used by us to re-analyze the H2 
absorption in the MS, we measure a S/N of ^ 6 per 
resolution element at 1020 A. Only for the wavelength 
range 1050 — 1082 A does the S/N ratio rise up to a 




Galactic longitude 

Fig. 1. — Hi column- density map of the Magellanic Stream in the 
direction of Fairall 9, based on 21cm data from GASS. The map 
shows the distribution of neutral gas in the LSR velocity range be- 
tween -j-lOO and -1-250 kms~^. The position of Fairall 9 is indicated 
with the white ellipse, whose shape reflects the beam size of the 
radio observations from GASS. 

maxium of ^ 15 per resolution element due to the in- 
creased background flux from the broad, redshiftcd Ly/3 
emission from Fairall 9. 

2.4. GASS 21cm observations 

The Hi 21cm data for the Fairall sightline were taken 
from the Galactic-All Sky Survey (GASS, McClure- 
Grifliths et. al.2009, Kalberla et al.2010). The survey 
was observed with the 64-m radio telescope at Parkes. 
The data cubes have an angular resolution of 15. '6, lead- 
ing to an RMS of 57 mK per spectral channel (Av — 
0.8 kms"^). This value translates to an Hi column den- 
sity detection hmit of A^(Hi)iini = 4.1 x 10^® cm~^, as- 
suming a Gaussian-like emission line with a width of 
20 km s^ FWHM. Fig. 1 shows a H i column density map 
of the local environment of the Magellanic Stream cen- 
tered on Fairall 9. 

2.5. Spectral analysis methods 

Our strategy for the analysis of the optical and UV 
absorption-line data of Fairall 9 combines different tech- 
niques to optimally account for the different spectral res- 
olutions and S/N in the data. The reduced and coad- 
ded spectra from UVES, COS, and FUSE first were 
continuum-normalized using low-order polynomials that 
were were fitted locally to the spectral regions of interest. 

For the high-resolution VLT/UVES data we have used 
Voigt-profile fitting to decompose the MS absorption pat- 
tern in the optical lines of Can, Nai, and Tin into indi- 
vidual absorption components (Voigt components) and 
to derive column densities (N) and Doppler parameters 
(6 values) for the individual components. For the fitting, 
we made use of the FITLYMAN package implemented in 
the ESO-MIDAS analysis software (Fontana & Ballester 
1995). Laboratory wavelengths and oscillator strengths 
have been taken from the compilation of Morton (2003). 
Note that with our fitting technique we are able to mea- 
sure b values smaller than the intrumental resolution as 
we simultaneously fit the line doublets of Can and Nai, 
so that relative strengths of the lines are taken into ac- 
count for the determination of both b and N. From 
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Fig. 2. — Optical absorption profiles of Call, Nal, and Till from 
VLT/UVES data of Fairall 9 are shown. Absorption from the Mag- 
ellanic Stream is seen at LSR velocities between +130 and +240 
kms~^. The solid red line display the best possible Voigt profile fit 
to the data. The individual velocity components are indicated by 
the tic marks. In the upper panel, the Hi 21cm emission spectrum 
toward Fairall 9 from the GASS data is plotted for comparison. 

this fitting procedure we obtain a component model for 
Fairall 9 sightline, in which the LSR velocity centroids of 
the individual absorption components and the b values 
for the low ions are defined. 

In addition to Voigt-profile fitting, we have used the 
apparent optical depth method (AOD method; Savage 
& Sembach 1991) to derive total gas column densities 
for the (unsaturated) optical absorption profiles of Can, 
Nai, and Tin. The AOD analysis was made using the 
custom-written MIDAS code span that allows us to mea- 
sure equivalent widths and AOD column densities (and 
their errors) in absorption spectra from a direct pixel 
integration. 

For the medium-resolution HST/COS data we have 
used profile modeling and the AOD method to derive 
column densities and column density limits for the 
various different low, intermediate, and high ions that 
have detectable transitions in the COS wavelength 
range. The ion transitions in the COS data consid- 
ered in this study include C ii A1334.5, Cii"^ A1335.7, 
Civ AA1548. 2, 1550.8, Ni AA1199.6, 1200.2, 1200.7, 
Nv AA1238. 8, 1242.8, Oi A1302.2, Alii A1670.8, 
Sill AA1190.4, 1193.3, 1260.4, 1304.4,1526.7, Sim 
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Fig. 3. — Velocity profiles of various low ions in the COS data 
of Fairall 9 are shown. The FUV data for Sill A1020.7 are from 
FUSE. The velocity components related to gas in the Magellanic 
Stream, as identified in the optical Call data, are indicated with the 
tic marks. The red solid line indicates the best-fitting absorption 
model for the UV absorption, as described in Sect. 3.2. The single, 
saturated lines of Cll A1334.5 and Ol A1302.2 are not considered 
in our model. 

A1206.5, Si IV AA1393.8, 1402.8, Pii A1152.8, Sii 
AA1250.6, 1253.8, 1259.5, Fell AA1143.2, 1144.9, 1608.5, 
and Niii A1370.1. For the profile modeling of neutral 
and singly-ionized species that trace predominantly 
neutral gas in the MS we have used as input the compo- 
nent model defined by the optical Can absorption (see 
above). In this model, the LSR velocities and b values 
of the seven absorption components are constrained by 
the Can fitting results, while the column density for 
each component is the main free parameter that can be 
varied for each ion listed above. Our previous studies 
of Can in the Galactic halo (Richter et al. 2005, 2009; 
Ben Bekhti et al.2008, 2011; Wakker et al. 2007, 2008) 
and in intervening absorbers at low redshift (Richter et 
al. 2011) have demonstrated that Can is an an excellent 
tracer for the distribution of neutral and partly ionized 
gas and its velocity-component structure, even in regions 
where Can is not the dominant ionization state. 

Based on the Can model, and using a modified ver- 
sion of FITLYMAN, we have calculated for each individual 
UV line a synthetic absorption profile, for which we have 
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Fig. 4. — COS velocity profiles of the intermediate and high ions 
Sim, Civ and Silv. The velocity structure in these ions is very 
different from that of the low ions, suggesting that they trace spa- 
tially distinct regions within the MS. High-velocity Si IV absorption 
between -f 150 and -f200 kms~l in the Si IV A1402.8 Une is possibly 
blended with intergalactic absorption. 

convolved an initial Voigt profile with input parameters 
{vi,bi,Ni) for i = 1...7 with the COS line-spread func- 
tion (LSF) that is appropriate for the wavelength of the 
line. For the COS LSF we have used the improved LSF 
model described by Kriss (2011). Then, the columns 
Ni were varied to minimize the differences between the 
synthetic absorption profile and the observed COS data. 
This method delivers reliable (total) column densities for 
those ions that have multiple transitions with substan- 
tially different oscillator strengths in the COS wavelength 
range and for individual lines that are not fully saturated. 
More details about the accuracy of this method are pre- 
sented in Sect. 3.1 and in the Appendix. 

For the intermediate and high ions (Si ill, Si iv, C iv) 
the component structure (and b values) are expected to 
be different from that for the low ions (as the gas phase 
traced by these ions often is spatially distinct from the 
gas phase traced by the low ions), but no information 
is available on the true component structure of these 
ions from the optical data. Therefore, we did not try to 
model the absorption profiles of the high-ion lines, but 
estimated total column densities (and limits) for these 
ions solely from the AOD method. Similarly, we used 
solely the AOD method to determine the column density 
of Cii*. 

For the analysis of H2 in the Magellanic Stream de- 
tected in the FairaU9 FUSE data (Richter et al. 2001a), 
we have modeled the H2 absorption using synthetic spec- 
tra generated with FITLYMAN. As model input we take 



into accout the component structure and line widths 
seen in the optical Call/Nal absorption, together with 
a Gaussian LSF according to the FUSE spectral resolu- 
tion. H2 wavelengths and oscillator strengths have been 
adopted from Abgrall & Roueff (1989). Details on the 
H2 modeling are presented in Sect. 3.3. 

3. COLUMN-DENSITY MEASUREMENTS 
3.1. Metal absorption in the optical 

Optical absorption related to gas in the Magellanic 
Stream in the velocity range 120 — 230 kms""'^ (this veloc- 
ity range is defined by the H i 21 cm data of the Stream) 
is detected in the UVES spectrum of Fairall 9 in the 
lines of Can (AA3934.8, 3969.6), Nai (AA5891.6, 5897.6), 
and Till (A3384.7; 3Aa detection). Fig. 2 shows the 
normalized absorption profiles of these ions plotted on 
the LSR velocity scale together with the Hi 21cm emis- 
sion profile from CASS. Can absorption is detected 
in seven individual absorption components centered at 
^^LSR = +143, +163, +172, -1-188, -1-194, -1-204 and +218 
kms""'^ with column densities in the range log N{Ca,u) = 
11.29 — 11.68. The red solid line in Fig. 2 indicates the 
best Voigt-profile fit to the data. All Can, Nai and Tin 
column-density measurements are summarized in Table 
1. 

The strongest Can absorption component in the MS 
(in terms of the absorption depth) is component 4 at 
+ 188 kms""'^; this very narrow component is also de- 
tected in both Nai lines (Fig. 2). From the simultaneous 
fit of the Ca n and Na i doublets we obtain a very small 
Doppler parameter for component 4 of 6 = 1.8 kms^^. 
The small b value for this component is confirmed by fits 
to the individual lines of Can and Nai, which all imply 
6 < 2 kms~^. The detection of Nai together with the 
small b value indicates that the gas in component 4 is 
relatively cold and dense and possibly is confined in a 
dense core with little turbulence. Components 3,4, and 
5 are also detected in Tin (at SAa significance; Fig. 2). 
Since Ti n and H i have almost identical ionization po- 
tentials (Table 2), the detection of Tin suggests that 
most of the neutral gas column density is contained in 
these three components. This conclusion is supported 
by the Hi 21cm emission profile from CASS (Fig. 2; up- 
per panel), which also has its maximum in the velocity 
range between +170 and +200 kms~^. Weak and rel- 
atively broad (6 — 7.4 kms^^) Nai absorption is also 
detected in component 1 at +143 kms~^. 

As mentioned in Sect. 2.5, we also have used the AOD 
method to determine the total column densities for Can, 
Nai and Tin by integrating over the velocity range rel- 
evant for MS absorption (i.e., +130 kms^^ < i^lsr < 
+230 kms~^.) The column densities derived by these 
two different methods agree very well within their la 
error ranges (see Table 1, rows 10 and 11). 

3.2. Metal absorption in the UV 

In Fig. 3 we show normalized UV absorption profiles 
of the low ions Cii, Ni, Oi, Aln, Sin, Sii, and Fen. 
The absorption profiles of the intermediate and high ions 
Sim, Civ and Siiv are shown in Fig. 4. Following the 
procedure described in Sect. 2.5, we have reconstructed 
the absorption pattern of UV metal lines of the low ions, 
using the component model defined by the optical lines of 
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TABLE 1 

Summary of ion column-density measurements 





Component 


1 


2 


3 


4 


5 


6 


7 








V [kms^-'-] 


+ 143 


+163 


+ 172 


+188 


+194 


+204 


+218 








b [kms-i] 


7.4 


7.9 


6.4 


1.8 


3.4 


4.9 


5.1 








o-f, [kms^i] 


2.4 


4.1 


2.2 


0.3 


1.4 


1.7 


2.0 






Ion 


Instrument 


logAfi 


logAr2 


logNs 


logAr4 


logAfs 


logA^e 


logA^T 


logA^tot.flt 


logAftot.AOD 


Hi 


CASS 


















19.95 


Nal 


VLT/UVES 


10.92 






11.05 








11.29 + 0.04 


11.33 + 0.05 


Can 


VLT/UVES 


11.58 


11.29 


11.63 


11.68 


11.58 


11.43 


11.35 


12.37 + 0.02 


12.38 + 0.03 


Tin 


VLT/UVES 






11.11 


11.15 


11.24 






11.27 + 0.06 


< 12.00 


Ci 


HST/COS 


















< 13.02 


Cii 


HST/COS 


















> 14.93 


Cii* 


HST/COS 


















13.35 + 0.07 


Civ 


HST/COS 


















13.73 + 0.04 


Ni 


HST/COS 


13.65 


13.42 


13.84 


14.10 


13.98 


13.54 


13.46 


14.63 + 0.05 


> 14.51 


Oi 


HST/COS 


















> 14.77 


Alii 


HST/COS 


12.55 


12.53 


12.92 


12.58 


12.69 


12.37 


12.24 


13.45 + 0.06 


> 13.18 


Sin 


HST/COS 


14.10 


13.75 


14.17 


14.13 


14.24 


13.92 


13.89 


14.90 + 0.04 


15.01 + 0.17 


Sim 


HST/COS 


















> 13.57 


Si IV 


HST/COS 


















13.03 + 0.06 


Pii 


HST/COS 


















< 12.87 


Sii 


HST/COS 


13.95 


13.70 


14.03 


14.00 


14.09 


13.77 


13.74 


14.77 + 0.02 


14.77 + 0.04 


Fell 


HST/COS 


13.47 


13.61 


13.94 


13.85 


13.94 


13.67 


13.65 


14.61 + 0.04 


14.59 + 0.04 


Niii 


HST/COS 


















< 13.37 



Can. As already mentioned, the modeling method pro- 
vides relevant results only for those ions that have mul- 
tiple transitions in the available COS vi^avelength band 
and for single lines that are not (or at most mildly) sat- 
urated. In our case, the modeling method could be used 
to determine column densities for Ni, Alii, Sin, Sii, 
and Fell, while for the fully saturated (single) lines of 
C II and O i the modeling does not yield relevant column 
density limits. In Fig. 3, the best-fitting column-density 
model for each ion is overlaid with the red solid line. The 
best-fitting model assumes the same b values for the indi- 
vidual subcomponents as derived for Ca ii from the Voigt 
fit of the UVES data. The ion column densities for each 
velocity component are listed in Table 1 in rows 2 — 7; 
the total column density for each ion is given in row 10. 
As can be seen in Fig. 3, the absorber model successfully 
reproduces the shape of the absorption lines in the UV. 
From Table 1 follows that the relative column-densities 
in the seven velocity components differ slightly from ion 
to ion. These differences are expected, because most ions 
have different ionization potentials and therefore trace 
slightly different gas phases within the absorber. In ad- 
dition, differential dust depletion for some elements may 
also be relevant in this context. Because of the relatively 
limited spectral resolution of the COS instrument and 
the resulting overlap between neighbouring velocity com- 
ponents the accuracy for the determination of the column 
densities for individual subcomponents is relatively low 
(~ 0.30 dex, typically). Much better constrained are the 
total column densities (Table 1, row 10) for the ion multi- 
plets Ni, Sii, Sin and Fcii, which have Icr uncertainties 
< 0.05 dex. This is because the integrated (total) col- 
umn density for each of these ions is determined predom- 
inantly by the integrated equivalent widths of the weaker 
lines and the equivalent width ratios of the weaker and 
stronger transitions. In the Appendix we present addi- 
tonal model plots to provide a more detailed insight into 



the allowed parameter range (bi,Ni) for the individual 
subcomponents and their impact on the shape of the ab- 
sorption profiles and total column-density estimates. 

Next to the absorption modeling wc have used the 
AOD method to determine total column densities (or 
column-density limits) in the COS data for each of the 
ions listed above (Table 1, row 11). For the unsaturated 
lines of Sii, Sin, and Fell the total column densities de- 
rived from the AOD method arc in excellent agreement 
(< Icr) with the total column densities determined from 
the profile modeling. Note that the absorption of C ii* 
in the MS will be discussed separately in Sect. 4.4. 

3.3. Hi 21cm emission 

The 21cm velocity profile shown in the upper panel 
of Fig. 2 indicates H i emission from gas in the Magel- 
lanic Stream in the LSR velocity range between -1-100 
and -1-250 km~^. The MS emission shows an asym- 
metric pattern with a peak near wlsr ~ +180 km~^, 
surrounded by weaker satellite componets at lower and 
higher radial velocities. The overall shape of the emission 
profile mimics that of the inversed absorption profiles 
of the low-ion lines of N i and S ii indicating that both 
the pencil-beam absorption data and the CASS emission 
data, that have a beam-size of 15.'6, sample the same 
physical structure. Integrating the GASS data in the 
above given velocity range for the single pixel that cov- 
ers the position of Fairall 9 we obain a total H i column 
density in the Stream of log TV (Hi) = 19.95 (Nijii) = 
(8.96 ±0.09) X 10^^ cm"^). The pixel value is somewhat 
smaller than the beam-averaged column density, which is 
log AT (Hi) = 19.98. For comparison, Gibson et al. (2000) 
find log 7V(Hi) = 19.97 (Ar(Hi= (9.35 ± 0.47) x lO^^ 
cm~^) using Parkes 21cm observation of Fairall 9 with 
a beam size of 14'0. These value agree with each other 
within their la error ranges. From the 21cm data of 
the Leiden- Argentine- Bonn Survey (LAB; Kalberla et 
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al. 2005) with ^ 0.5 deg resolution we obtain for the same 
velocity range an H i column density in the Stream of log 
7V(Hi) = 19.97, thus very similar. 

For the determination of metal abundances (using the 
GASS data in combination with the pencil-beam UVES 
and COS data) we adopt the value of log Af(Hi) = 19.95 
for the total H i column density in the Stream toward 
Fairall 9, together with an appropriate estimate of the 
systematic error of 0.03 dex that accounts for the lim- 
ited spatial resolution of the 21cm data. As it is shown 
in the Appendix, this H i column density in the MS is 
in agreement with with the observed shape of the Ly a 
absorption line in the COS data of Fairall 9, which pro- 
vides an independent (albeit less stringent) measure for 
the neutral-gas column density in the Stream in this di- 
rection (see also Wakkcr, Lockman, & Brown 2011 for a 
detailed discussion on this topic). 

3.4. Molecular hydrogen absorption 

The detection of H2 absorption in the Magellanic 
Stream in the FUSE data of Fairall 9 has been reported 
previously by Richter et al. (2001a) and Wakker (2006). 
H2 absorption is seen near wlsr = +190 kms~^ in the 
rotational levels J = 0, 1 and 2. Although the overall 
S/N in the FUSE data is low (< 15 per resolution ele- 
ment; see Sect. 2.5), the H2 absorption is clearly visible in 
many lines due to the relatively high H2 column density 
(Richter et al. 2001a). 

We have reanalyzed the Fairall 9 FUSE data based on 
an improved data reduction pipeline and the component 
model discussed in the previous sections. While Richter 
et al. (2001) have used a curve-of-growth technique to 
derive the column densities N(J) and b values of the H2 
rotational states J = 0, 1, 2, assuming as single MS ab- 
sorption component, we here model the absorption spec- 
trum of the MS H2 absorption in the FUSE FairaU9 
data using synthetic spectra. For our model we assume 
that the H2 absorption arises in the region with the high- 
est gas density in the Stream, which is component 4 at 
I'LSR = +188 kms~^ detected in the UVES data in Can, 
Nal, and Tin (Table 2). In fact, from absorption studies 
of Nai and H2 in interstellar gas it is well known that 
these two species trace the same gas phase (i.e., cold 
neutral and cold molecular gas; e.g., Welty et al. 2006). 

We adopt 5turb = 1-8 km~^ for the H2 absorption in 
component 4, as suggested by the Nal and Can absorp- 
tion in the same component (Table 2). We further as- 
sume (and prove later) that the kinetic temperature of 
the H2 absorbing gas is Tkin ~ 100 K in component 4. It 
is usually assumed that the Doppler parameter is com- 
posed of a turbulent and a thermal component, so that 
6^ = 6j^j.jj + 6j|j. For light elements and molecules (such 
as H2) 5th must not be neglected. We therefore adopt 
^(H2) = 2.0kms~^ for our H2 model spectrum. 

Fig. 4 shows our best fitting H2 model (red solid line) 
together with the FUSE data for the wavelength range 
between 1070 and 1080 A. Since the 6 value is fixed to 
^(H2) — 2.0kms~^ and the velocity to wlsr = +188 
kms~^, the only free parameters in our model are the 
II2 column densities N{J). Table 2 summarizes the best 
fitting II2 column densities using this method. The 
given errors represent the uncertainties for log N{J) un- 
der the assumption that b = 2.0 kms~^. The result- 



TABLE 2 
Summary of H2 column-density modeling 



b^ [kms ^] 



log N{J) 



Tj [K]b^ 






2.0 


17.53 ±0.24 


^•^-39 


1 


2.0 


17.68 ±0.26 


qQ+149 


2 


2.0 


16.46 ±0.21 


<172 


3 


2.0 


< 15.60 


< 172 


Sum 




ivqq+O.lS 





^ fe- value fixed to 2.0 kms ^; see Sect. 3.4. 

^ Equivalent Boltzmann temperature implied by N{J); see 

Sect. 4.4.2. 

ing total H2 column density is log iV(H2) = 17.93^0.16 
and the fraction of hydrogen in molecular form is log 
/h, =log[27V(H2)/(iV(Hi) +2Af(H2))] = -2.03. Note 
that the value for log /hs derived in this way represents 
a sightline-average; the local molecular fraction in the II2 
absorbing region must be higher (because A^(Hi) in this 
region is smaller). 

The II2 column densities we derive from our spectral 
modeling are substantially (~ 1 dex) higher than the val- 
ues presented in the earlier studies of the Fairall 9 FUSE 
data (Richter et al. 2001a; Wakker 2006). The reason for 
this discrepancy is the substantially lower b value that we 
adopt here. This more realistic b value for the H2 absorb- 
ing gas is a direct consequence of the previously unknown 
complex absorption pattern in the MS toward Fairall 9, 
that is visible only in the high-resolution optical data of 
Can and Nai (see Sect. 3.1). Our study suggests, that, 
without the support of high-resolution and high S/N op- 
tical or UV data, that can be used to resolve the true 
sub-component structure of the absorbers and to deliver 
reliable b values for the subcomponents, the interpreta- 
tion of interstellar II2 absorption lines in low-resolution 
and low-S/N data from FUSE can be afflicted with large 
systematic uncertainties. 

4. PHYSICAL AND CHEMICAL CONDITIONS IN THE GAS 
4.1. Ionization conditions and overall metallicity 

The presence of high, intermediate and low ions, as 
well as molecular hydrogen, in the Magellanic Stream 
toward Fairall 9 in seven individual velocity components 
implies a complex multi-phase nature of the gas. To 
obtain meaningful results for the chemical composition 
and physical conditions from the measured ion column 
densities, it is necessary to consider in detail the ion- 
ization conditions in the absorbing gas structure. For 
an element M we define its relative gas-phase abundance 
compared to the solar abundance in the standard way 
[M/H] == log (M/H) - log (M/H)©. Solar reference abun- 
dances have been adopted from Asplund et al. (2009). 

The best ion to study the overall metallicity of neu- 
tral and weakly ionized interstellar gas is O i. Neutral 
oxygen and neutral hydrogen have similar ionization po- 
tentials and there is a strong charge-exchange reaction 
that couples both ions in gas with sufficiently high den- 
sity. For our sightline, however, the O i column density is 
not well constrained because of the saturation of the O i 
A1302.2 line and the lack of other, weaker O I transitions 
in the COS wavelength range. Because it is impossible 
to measure the O i column density in the MS using our 
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Fig. 5. — Continuum-normalized FUSE spectrum of Fairall9 in the wavelength range between 1075 and 1081 A. Absorption by H2 from 
the rotational levels J = 0, 1 and 2 related to gas in the Magellanic Stream are indicated above the spectrum using the common notation 
scheme. The best-fitting model for the H2 absorption is overplotted with the red solid line. 



component models and the AOD method to a satisfac- 
tory accuracy, we do not further consider the Oi A1302.2 
line for our analysis. 

For the determination of the overall metallicity of the 
gas in the Magellanic Stream toward Fairall 9 we choose 
instead S 11, which is another useful ion for measuring in- 
terstellar abundances, as singly-ionized sulfur is an excel- 
lent tracer for neutral hydrogen without being depleted 
into dust grains (e.g.. Savage & Sembach 1996). We also 
have a very accurate column-density determination of 
log iV(Sll) = 14.77 ± 0.02 in the MageUanic Stream, 
based on the absorption modeling of the two Sii lines 
at 1250.6 and 1253.8 A and (independently) from the 
AOD method (see Table 1). Our value is slightly higher 
than the S 11 column density of log N{S 11) = 14.69 ± 0.03 
derived by Gibson et al. (2000) from HST/GRRS data. 

From the measured (Sii/Hi) ratio we obtain an initial 
estimate of the sulfur abundance in the Stream toward 
Fairall 9 of [S/H] « —0.30. Since the ionization potential 
of Sii (23.2 eV; see Table 3) is higher than that of Hi 
(13.6 eV), we have to model the ionization conditions in 
the gas to obtain a more precise estimate for the sulfur 
abundance (and the abundances of the other detected 
metals) in the Stream toward Fairall 9. While previous 
studies of Galactic halo clouds (Wakker et al. 1999; Gib- 
son et al. 2000) have demonstrated that the (S ii/H i) ra- 
tio is a very robust measure for the sulfur-to-hydrogen 
ratio, hardly being affected by the local ionizing radi- 
ation field, the ionization modeling is important to ob- 
tain a relibale estimate of the systematic uncertainty that 
arises from the use of S ii as reference indicator for the 
overall a abundance in the gas. For this task, we have 
used the photoionization code Cloudy (vlO.OO; Ferland et 
al. 1998), which calculates the expected column densities 
for different ions as a function the ionization parameter 
U — riy/nii for a gas slab with a given neutral gas col- 
umn density and metallicity, assuming that the slab is 
illuminated by an external radiation field. As radiation 
field we adopt a combined Milky Way plus extragalactic 
ionizing radiation field based on the models by Fox et 
al. (2005) and Bland-Hawthorn & Maloney (1999,2002), 
appropriately adjusted to the position of the MS relative 
to the Milky Way disk. 

While our absorption modeling provides rough esti- 
mates for the ion column densities in the seven sub- 



components, we refrain from trying to model the ion- 
izations conditions in the individual components. This 
is because the unknown Hi column densities in these 
components, together with the unknown geometry of 
the overall gas structure, would lead to large systematic 
uncertainties for such models. Instead, we have used 
Cloudy to obtain integrated elemental abundances in the 
gas of the Magellanic Stream toward Fairall 9 (represent- 
ing the optical-depth weighted mean of the individual 
element abundances in all subcomponents). It is, how- 
ever, important to emphasize that meaningful results for 
integrated elemental abundances from Cloudy can be ob- 
tained only for those ions for which the column densities 
in the ionization model do not strongly depend on the 
gas density (and ionization parameter), because the lat- 
ter quantities are expected to vary substantially among 
the individual subcomponents. 

For our Cloudy model we assume a neutral gas col- 
umn density of log A^(Hi) = 19.95, an overall metallicity 
of [Z/H] = —0.3, and solar relative abundances of the 
metals from Asplund et al. (2009), based on our results 
described above. In Fig. 5, upper panel, we display the 
expected logarithmic column densities for low and inter- 
mediate ions as a function of log U and log nn for this 
model. The measured total column densities for these 
ions (Table 1, 10th column) are plotted on the right hand 
side of the panel. The total Ca ii column density of log 
A^(Caii) = 12.37 sets a limit for the (averaged) ionisation 
parameter and density in the gas of log (U) < —3.55 and 
log (nn) > —1.95. The corresponding density range that 
we consider as relevant for our Cloudy model (2.5 <log 
TT-H < —1.95) is indicated in Fig. 5 with the gray-shaded 
area. For this range, the expected column densities of 
Ni, Sin, Fell, Sli, Alii, Niii, and Tin are nearly inde- 
pendent of log U and log nn- For Sii, in particular, the 
column density varies within only 0.03 dex in the above 
given density range, supporting the previous conclusions 
from Wakker et al. (1999) and Gibson et al. (2000). In the 
lower panel of Fig. 5 we show the differences (A7V) be- 
tween the measured column densities and the mean col- 
umn densities predicted by Cloudy. For Niii and Can, 
the arrows indicate the upper limit and the allowed range 
for AA^, respectively. 

From the comparison between the predicted and mea- 
sured ion column densities, together with the solar ref- 
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Fig. 6. — Upper panel: Cloudy photoionization model of low- 
ionization species in the MS toward Fairall 9. The colored lines 
indicate the predicted ion column densities as a function of the 
ionization parameter U and gas density njj for a gas slab with a 
(total) neutral hydrogen column density of log Af(Hl) = 19.95 and 
a metallicity of [Z/H] = —0.30. Observed ion column densities are 
indicated on the right-hand side with the filled colored circles. The 
gray-shaded area marks the range for U and njj that is relevant for 
our study based on the oberved Call column density (see Sect. 4.1). 
Lower panel: differences between the observed and the predicted 
ion column densities, AN = A^cioudy ~ ^obs.i these differences are 
discussed in Sect. 4.1. 

erence abundances from Asplund et al. (2009), we obtain 
the following gas-phase abundances: [N/H] = — 1.15 ± 
0.06, [Si/H] = -0.57 ± 0.06, [Fe/H] = -0.86 ± 0.06, 
[S/H] = -0.30 ± 0.04, [Al/H] = -0.92 ± 0.08, [Ti/H] = 
-1.62 ± 0.07, [Ni/H] < -0.81, and [Ca/H] < 0. The 
listed errors include the uncertainties from the column 
density measurements for the metal ions and H i (Table 
1), the systematic uncertainty due to the beam size of 
the 21cm measurements (Sect. 3.3), and the systematic 
uncertainty for the ionization correction from the Cloudy 
model (see above). We here do not consider the system- 
atic errors that come with the solar reference abundances 
(Asplund et al. 2009; see Table 3, third row), because 
these errors are irrelevant for the comparison between 
our results and other abundance measurements, if the 
identical reference values are used. All gas-phase abun- 
dances derived with Cloudy are listed in Table 3, sixth 
row. 

The measured S/H ratio of [S/H] = -0.30 ±0.04 corre- 
sponds to an overall sulfur abundance in the Magellanic 
Stream toward Fairall 9 of 0.50^0 04 solar. The value for 
[S/H] is 0.25 dex higher than the value derived for the MS 
toward Fairall 9 by Gibson et al. (2000) based on GHRS 
data. Several factors contribute to this discrepancy: (a) 



the somewhat higher (-1-0.08 dex) Sii column density 
that we derive from the COS data, (b) the slightly lower 
Hi column density (—0.02 dex) that we obtain from the 
CASS data, and (c) the substantially higher (0.15 dex) 
solar reference abundance of sulfur (Asplund et al. 2009) 
that we use to derive [S/H] compared to the value from 
Anders & Grevesse (1989) used by Gibson et al. (2000). 

Remarkably, the sulfur abundance we derive here is 
higher than the present-day stellar sulfur abundances 
in the SMC ([S/H] = -0.53 ± 0.15) and in the LMC 
([S/H] = -0.42 ± 0.09; RusseU & Dopita 1992), but 
matches the present-day interstellar S abundances found 
in Hii regions in the Magellanic Clouds (Russell & Do- 
pita 1990). This interesting result will be further dis- 
cussed in Sect. 6. Because of the large abundance scatter 
of sulfur in the Magellanic Clouds and because S ap- 
pears to be underabundant compared to oxygen in the 
LMC and overabundant compared to oxygen in the SMC 
(Russel & Dopita 1992), the measured (S/H) ratio does 
not provide a constraining parameter to pinpoint the ori- 
gin of the gaseous material in the Stream toward Fairall 9 
in cither SMC or LMC. 

From our Cloudy model it follows that all other ele- 
ments listed above have gas-phase abundances (or abun- 
dance limits) that are lower than that of S. Next to the 
intrinsic nucleosynthetic MS abundance pattern, dust de- 
pletion is another important effect that contributes to 
deficiency of certain elements (e.g., Ti, Al, Fe, Si) in 
the gas phase. This aspect will be discussed in Sect. 4.3. 
Note, again, that we do not attempt to constrain a single 
value for log U or log nn with our Cloudy model, as our 
sightline passes a complex multiphase structure that is 
expected to span a substantial range in gas densities an 
ionization parameters. 

In this paper, we do not further analyze the gas phase 
that is traced by the intermediate and high ions Sim, 
C IV and Si iv, which probably is tracing the ionized 
boundary layer between the neutral gas body of the 
Stream and the ambient hot coronal halo gas of the Milky 
Way (see Fox et al. 2010). For completeness, we show in 
Fig. 4 the velocity profiles of the available transitions of 
Sim, Civ and Siiv in the COS spectrum of FairaU9. 
The velocity structure in these ions clearly is different 
from that of the low ions, suggesting that they trace spa- 
tially distinct regions in the MS. The high- ion absorption 
toward Fairall 9 will be included in a forthcoming paper 
that will discuss those ions in several lines of sight inter- 
secting the Magellanic Stream (A.J. Fox et al. 2014, in 
preparation) . 

4.2. Nitrogen abundance 

The measured gas phase abundance of nitrogen in the 
Magellanic Stream toward Fairall 9 is [N/H] — — 1.15 ± 
0.06, which is 0.85 dex lower than the abundance of sul- 
fur. In diffuse neutral gas the depletion of nitrogen into 
dust grains is expected to be very small (Savage & Sem- 
bach 1996). Consequently, the low nitrogen abundance 
in the Stream reflects the nucleosynthetic enrichment his- 
tory of the gas and thus provides important information 
about the origin of the gas. 

For the SMC, Russeh & Dopita (1992) derive a 
mean present-day stellar nitrogen abundance of [N/H] = 
— 1.20. This value is close to the value derived for the 
MS toward Fairall 9, although the gas in Stream has not 
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TABLE 3 

Summary of integrated gas-phase abundances and dust depletion values 



Ion (X) 


Z 


LP.'' 

[eV] 


log(M/H)b 
+12 


log(X/Hi) - 


[M/H]d 


log 5(M) 


Cl 


6 


11.26 


8.43 ±0.05 


< -6.93 






Cii 


6 


24.38 


8.43 ± 0.05 


> -5.02 






Ni 


7 


14.53 


7.83 ± 0.05 


-5.32 ±0.05 


-1.15 ±0.06 




Oi 


8 


13.62 


8.69 ± 0.05 


> -1.87 






Nal 


11 


5.14 


6.24 ± 0.04 


-8.66 ±0.01 






Alii 


13 


18.83 


6.45 ± 0.03 


-6.50 ±0.06 


-0.92 ±0.08 


-0.62 


Sin 


14 


16.35 


7.51 ± 0.03 


-5.05 ±0.04 


-0.57 ±0.06 


-0.27 


Pii 


15 


19.77 


5.41 ± 0.04 


< -7.08 






Sii 


16 


23.34 


7.12 ± 0.03 


-5.18 ±0.02 


-0.30 ±0.04 


(0) 


Can 


20 


11.87 


6.34 ± 0.04 


-7.58 ±0.02 


< -0.30 


to -1.58 


Tin 


22 


13.58 


4.95 ± 0.05 


-8.68 ±0.06 


-1.63 ±0.07 


-1.33 


Fell 


26 


16.12 


7.50 ± 0.04 


-5.34 ±0.04 


-0.86 ±0.06 


-0.56 


Niii 


28 


18.17 


6.22 ± 0.04 


< -6.58 


< -0.81 


< -0.51 



^ I. P.= ionization potential. 

^ Solar reference abundance for element M from Asplund et al. (2009). 

'^ Listed errors include uncertainties from the column-density measurements of X and Hi. 

•^ Gas phase abundance for element M, defined as [M/H] = log (M/H) - log (M/H)© , derived from the CLOUDY model. The listed 

errors include uncertainties from the column-density measurements, the beam-smearing of the 21cm data, and the ionization 

correction from the Cloudy modeling. 

<= Depletion value, defined as log S(M)= [M/H]ms-[S/H]ms- 



been enriched with nitrogen since it was stripped off its 
parent galaxy. The nitrogen abundances found in SMC 
H II regions and supernova remnants span a rather wide 
rangcof-1.37<[N/H] < -0.95 (Russell & Dopita 1990). 
Thus, the observed nitrogen abundance in the MS toward 
Fairall9 would be in accordance with an SMC origin of 
the gas only, if the (mean) nitrogen abundance in the 
SMC has not increased substantially after the Stream 
was separated ^1 — 2 Gyr ago. 

In the LMC, in contrast, the mean present-day stellar 
nitrogen abundance is found to be much higher ([N/H] — 
—0.69). The N abundance range in LMC H ii regions and 
supernova remnants is [N/H] = —0.86 to —0.38, the min- 
imum [N/H] still being substantially higher than the ni- 
trogen abundance in the Stream toward Fairall9. There- 
fore, an LMC origin of the MS is plausible only if the 
LMC has substantially increased its nitrogen abundance 
during the last ^ 2 Gyr. 

In Sect. 5.1 we will further consider these interesting 
aspects and discuss the observed N abundance in the MS 
in the context of the enrichment history of the Magellanic 
Clouds. 

4.3. Dust depletion pattern 

Heavy elements such as Al, Si, Fe, Ca, Ti, Fe and Ni 
are known to be strongly depleted into dust grains in in- 
terstellar gas in the Milky Way and other galaxies, while 
other elements sich as S, O, N are only very mildly or not 
at all affected by dust depletion (e.g.. Savage & Sembach 
1996; Welty et al. 1997). 

To study the depletion pattern of these elements in the 
Magellanic Stream we define the depletion values from 
the relative abundance of each depleted element, M, in 
relation to the abundance of sulfur. We define accord- 
ingly log (5(M)=[M/H]ms-[S/H]ms- Using this defini- 
tion, the depletions values are log S{Si) = —0.27, log 
(5(Fe) = -0.56, log S{Ti) ^ -1.32, log 6{Al) = -0.62, 
and log (5(Ni) < —0.51 (see also Table 3, seventh row). 
Such depletion values are typical for warm Galactic halo 



clouds, as derived from UV absorption-line measure- 
ments (see Savage & Sembach 1996; their Fig. 6). 

Note that the depletion values (or limits) would be 
different, if the relative chemical abundances of the de- 
pleted elements would differ from the relative chemical 
composition of the sun (Asplund et al. 2009). In view 
of the abundance patterns found in the LMC and SMC 
(Russell & Dopita 1990,1992), this is actually a likely 
scenario, but no further conclusions can be drawn at this 
point without knowing the exact intrinsic chemical com- 
position of the MS gas toward Fairall9. 

Since the expected Can column density in MS gas is 
expected to depend strongly on the local gas density 
and ionization parameter (Fig. 6), the depletion value 
of Ca cannot be tightly constrained, but may lie any- 
where in the range log S{Ca,) > —1.58. Under typi- 
cal interstellar conditions, Ca is strongly depleted into 
dust grains with large depletion values that are similar 
to those of other elements that have condensation tem- 
peratures above T = 1500 K (e.g., Ti; Savage & Sembach 
1996). If one assumes that log ^(Ca) «log S{Ti) for the 
MS toward Fairall9, then Fig. 6 would imply that the 
average density in the neutral gas in the Stream in this 
direction is log tih ~ or tih ~ 1 cm~^. This density, 
together with the neutral gas column density, would im- 
ply a characteristic thickness of the absorbing neutral gas 
layer of d = N(ai)/nii « 30 pc. 

At first glance, one may regard this as a relatively small 
value. The complex velocity-component structure of the 
absorber indicates, however, that the absorber is com- 
posed of several small cloudlets along the line of sight 
that (together) are extending over a spatial range that is 
larger than the value indicated by the mean thickness. In 
view of the relatively high metal abundance derived for 
the kinematically complex absorber toward Fairall9, the 
presence of a population of metal-enriched gas clumps at 
small scales in the MS is not an unlikely scenario. 
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Fig. 7. — COS velocity profiles of Cll* A1335.7 (upper panel). 
C II* absorption from the Magellanic Stream is visible between 
+120 and +200 kms~^. For comparison, the velocity profile of 
Sll A1250.6 is also shown (lower panel). The component structure 
of the low ions is indicated with the tic marks. 

4.4. Physical conditions in the H2 absorbing region 

As discussed in Sect. 3.1, the UVES spectrum of 
Fairall 9 shows the presence of Nal, Can, and possibly 
Tin in component 4. In the same component, the de- 
tected H2 absorption is expected to arise (Sect. 3.3), indi- 
cating the presence of a cold, dense (and predominantly 
neutral) gas phase in the MS in this direction. In the 
following, we want to combine the information from the 
different instruments to explore the physical and chemi- 
cal conditions in the II2 absorbing gas phase in the MS 
toward Fairall 9. 

4.4.1. Gas density 

While our Cloudy model does not provide any rele- 
vant information on the local gas densities in the individ- 
ual subcomponents, we can use the observed molecular 
hydrogen fraction in the gas to constrain nu in the II2 
absorbing component 4, assuming that the abundance of 
II2 relative to H i is governed by a formation-dissciation 
equilibrium. In a ormation-dissciation equilibrium the 
neutral to molecular hydrogen column density ratio in 
an interstellar gas cloud is given by 



A^(HI) _ (fc)/3 
N{R2) ~ i?n(HI). 



(1) 



where (fc) « 0.11 is the probability that the molecule 
is dissociated after photo absorption, (3 is the photo- 
absorption rate per second within the cloud, and R is 
the II2 formation rate on dust grains in units cm'^s^^. 
For low molecular hydrogen fractions we can write nn — 
n(III) + 2n(H2) w ri(HI). The parameter </> < 1 in equa- 
tion (1) describes the column-density fraction of the Hi 
that is physically related to the H2 absorbing gas, i.e., 
the fraction of the neutral gas atoms that can be trans- 
formed into H2 molecules (see Richter et al. 2003 for de- 
tails). Our absorption modeling of the undepleted low 
ions S II and N i indicates that component 4 contains 
^ 15 — 30 percent of neutral gas column, so that we 
assume = 0.15 — 0.30 in lack of any more precise infor- 
mation. 

For interstellar clouds that are optically thick in II2 
(i.e., for log iV(Il2) ^ 14) H2 line self-shielding needs to 
be considered in equation (1), because the self-shielding 
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Fig. 8. — Rotational excitation of H2 in the Magellanic Stream. 
The fit for the excitation temperatures (assuming a Boltzmann 
distribution) for J = 0, 1 and J = 2, 3 are indicated with the solid 
and dashed line, respectively. 



reduces the photo-absorption rate (/3) in the cloud inte- 
rior. Draine & Bertoldi (1996) find that the H2 self- 
shielding can be expressed by the relation /3 — 5/3o, 
where S = (A^H2/10"cm-2)-0-75 < 1 ig the self-shielding 
factor and $0 is the photo absorption rate at the edge of 
the cloud. The parameter /3o is directly related to the 
intensity of the ambient UV radiation field at the edge 
of the H2-bearing clump. Compared to the Milky Way 
disk, where UV bright stars dominate the mean photo- 
absorption rate of /3o,mw = 5.0 x 10~^° s~^ (e-g-, Spitzer 
1978), the value for /3o in the Magellanic Stream is ex- 
pected to be substantially smaller. This is because the 
solid angle of the Milky Way disk is relatively small at 
50 kpc distance and the contribution of the extragalactic 
UV background to /?o,ms is expected to be (relatively) 
small, too. From the model by Fox et al. (2005) follows 
that the UV flux between 900 and 1100 A is reduced by a 
factor of 160 compared to the mean flux within the Milky 



Way disk, so that we assume /3o,ms = 3.1 x 10 ^■^ s ^. 
With a total H2 column density of 8.5 x 10^^ cm~^ in the 
Stream the self-shielding factor S becomes 1.13 x 10""^, so 
that the photo-absorption rate in the cloud core is esti- 
mated as /3 = 3.5 X 10~^^ s~^ . For the H2 grain formation 
rate in the Magellanic Stream we adopt the value derived 
for the SMC based on FUSE H2 absorption-line data 
(Tumlinson et al.2002), i.e., Rms = 3 x 10"^*^ cm^s'^ 
This value is 10 times smaller than R in the disk of the 
Milky Way (Spitzer 1978). 

If we solve equation 3 for ^nn = (f>n(iil) and include 
the values given above, we obtain a density of (pnn w 1 
cm~^ or riH ~ 4 - 8 cm^^ for = 0.15 - 0.30. These 
densities are very close to the density derived for the 
II2 absorbing gas in the Leading Arm of the MS toward 
NGC3783 (riH « 10 cm'^; Sembach et al.2001). Note 
that II2 absorption has also been detected in the Magel- 
lanic Bridge (Lehner 2002). 

4.4.2. Gas temperature 
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In Fig. 8 we plot the measured logarithmie H2 eolunin 
density for each rotational state, log N{J), divided by 
its statistical weight, gj, against the rotational excitation 
energy, Ej. Rotational excitation energies and statistical 
weights have been adopted from the compilation of Mor- 
ton & Dinerstein (1976). The data points in Fig. 8 can be 
fitted by a Boltzmann distribution (i.e, a straight line in 
this plot), where the slope characterizes the (equivalent) 
excitation temperature for the rotational levels consid- 
ered for the fit. For the two rotational ground states 
J = 0, 1 the Boltzmann distribution thus is given by 
the equation N{1)/N{0) = gi/gQexp{-Eoi/kTQi). Us- 
ing this equation, together with our measured column 
densities for J = 0, 1, we find an excitation tempera- 
ture of Toi = 93^39^ K (Fig. 8 solid line). Because the 
two rotational ground states most likely are collisionally 
excited, Tqi represents a robust measure for the kinetic 
temperature in the H2 absorbing gas (e.g., Spitzer 1978); 
however, because of the saturation of the H2 lines and 
the resulting relatively large errors for iV(0, 1), the un- 
certainty for Tqi is substantial. In a similar manner, we 
derive for the rotational states J = 2,3 an upper limit 
for the equivalent Boltzmann temperature of r23 < 172K 
(Fig. 8 dashed line). This higher excitation temperature 
for J > 2 reflects excitation processes other than colli- 
sions, such as UV pumping and H2 formation pumping 
(Spitzer 1978). However, the relatively low temperature 
limit of r23 < 172 K reflects the low intensity of the am- 
bient dissociating UV radiation field at the position of 
the MS. 

The above given excitation temperatures deviate from 
our earlier estimates (Richter et al. 2001a). This is not 
surprising, however, since the values for N(J) have sub- 
stantially changed, too (see Sect. 3.3). The derived value 
of Tqi = 93 K lies in a temperature regime that is typi- 
cal for diffuse II2 gas in the disk of the Milky and in the 
Magellanic Clouds (e.g.. Savage et al. 1977; Tumlinson et 
al. 2002; de Boer et al. 1998; Richter et al. 1998; Richter 
2000). 

From the density and temperature estimate in the H2 
absorbing gas we now are able to provide a direct es- 
timate of the thermal pressure, P/k = nT, in the cold 
neutral gas in the MS. With rin = 8 cm^'' and T = 93 K 
we obtain a pressure of P/k = 744 cm^'^ K. This value is 
in good agreement with previous estimates for the ther- 
mal pressure in the cold neutral phase of the MS from 
Hi 21cm measurements (e.g., Wakker et al. 2002). 

The thickness of the H2 absorbing structure in the MS 
is d-a^ = (/)7V(Hi)/nH ~ 0.6 — 1.2 pc. This small dimen- 
sion of the absorbing gas clump explains the very small 
velocity dispersion (6 = 1.8 kms^^; see Table 1) that is 
measured for this component. 

4.5. Electron density and C^ cooling rate from Cll* 
absorption 

In Fig. 7 (upper panel) we show the LSR velocity pro- 
file of the Cii* A1335.7 line, based on the COS data. 
Weak absorption is visible at MS velocities, but the com- 
ponent structure is different from that of the other low 
ions (e.g., S ii; Fig. 7, lower panel). Because we could not 
identify any other (e.g., intergalactic) origin for this ab- 
sorption feature, we assume that it is caused by C 11* in 
the Magellanic Stream. The C 11* feature is characterized 



by a stronger and broader absorption peak near +150 
kms^^, falling together with components 1 and 2 defined 
in Table 1, and another weak and narrow component that 
is seen near +190 kms^^, most likely associated with 
component 4, in which H2, Nai, and Can is detected. 
Therefore, the relative strengths of the C 11* components 
appear to be inverted from those of the ground-state 
species. This indicates multi-phase gas, in which the ion- 
ization fractions and electron densities vary among the 
different absorption components. 

Under typical interstellar conditions, the relative pop- 
ulation of the fine-structure levels of ionized carbon (C''^) 
are governed by the balance between collisions with elec- 
trons and the radiative decay of the upper level into the 
ground state 2s'^2p^Pi/2. The Cii* A1335.7 transition 
arises from the 2s^2p^P-^/2 state, which has an energy of 
~ 8 X 10^ cV above the ground state. Measurements of 
the column-density ratios A^(Cii*)/7V(Cii) thus can be 
used to estimate the electron densities n^ in different in- 
terstellar environments, including Galactic high- velocity 
clouds (HVCs) and circumgalactic gas structures (e.g., 
Zech et al. 2008; Jenkins et al. 2005). 

For the C 11* absorption in the MS toward Fairall 9 we 
measure a total column density of log iV(C 11*) — 13.35 ± 
0.07 for the velocity range vlsr — 120 — 220 kms^^ using 
the AOD method. For the weak absorption associated 
with component 4 we derive log A^(C ii"^) = 12.63 ± 0.03 
(AOD). The C 11 1334.5 absorption in the MS is fully sat- 
urated and thus does not provide a constraining limit on 
iV(Cii) (see Table 1). We therefore use S 11 as a proxy, 
because the S ii/C 11 ratio is expected to be constant over 
a large density range (Fig. 6) and both elements are not 
expected to depleted into dust grains in the MS. Using 
a solar (S/C) ratio (Asplund et al. 2009), we estimate 
from the values listed in Table 1 that the total C 11 col- 
umn density in the MS is log A'tot(Cii) « 16.10 and the 
column density in component 4 is log A^4(Cii) « 15.30. 

Because the C+ fine-structure population depends 
strongly on the gas temperature (see, e.g., Spitzer 1978), 
we do not attempt to estimate a mean value for n(e) from 
A'^tot(C ll*)/iVtot(C 11), because we know that the temper- 
ature in this multi-phase absorber spans a large range be- 
tween ^ 100 K and probably few 1000 K, so that the re- 
sulting estimate for (rio) would be meaningless. Instead, 
we concentrate on component 4, for which we know the 
gas temperature and density from the analysis of the 
H2 rotational excitation (Sect. 4.4). Keenan et al. (1986) 
have calculated electron excitation rates for the C+ fine- 
structure transitions for a broad range of physical con- 
ditions in interstellar gas. Using our temperature and 
density estimates for the gas in component 4 (T = 93 K 
and riH =4 — 8 cm~^) together with their predicted pop- 
ulation rate ratios {2s^2p^P-ijil2s^2p^Pxji) for T = 100 
K and nn = 5 — 10 cm""^ (their Fig. 2), the measured 
column-density ratio log[iV(Cii*)/7V(Cii)] = -2.67 im- 
plies that the electron density in component 4 is small 
compared to nn, namely n^ < 0.05 cm~^. This low value 
for rio is in excellent agreement with expectations for a 
stable cold, neutral medium in gas with subsolar metal- 
hcities (Wolfire et al. 1995). 

Having a robust measure for A^tot(C 11*), it is also pos- 
sible to estimate the C+ cooling rate in the Magellanic 
Stream toward Fairall 9. Because ionized carbon is a ma- 



The Magellanic Stream toward Fairall 9 



13 



jor cooling agent of interstellar gas in a wide range of 
environments (e.g., Dalgarno & McCray 1972), the C+ 
cooling rate is an important quantity that governs the 
thermal state of diffuse gas inside and outside of galax- 
ies. For the hydrogen and electron densities in the Magel- 
lanic Stream (see above) the 0+ cooling rate is governed 
predominantly by spontaneous de-excitations, while col- 
lisional de-excitations can be neglected (e.g., Spitzer 
1978). Following Lehner, Wakker & Savage (2004), the 
C+ cooling rate per neutral hydrogen atom can be esti- 
mated as Ic = 2.89 X 10-2"iV(C+)/7V(Hi)ergs-i. If we 
adopt log N{Yi i) = 19.95 and log N{G ii*) = 13.35 we de- 
rive a mean (sightline average) C''" cooling rate per neu- 
tral hydrogen atom of log Ic = —26.14 for the MS toward 
Fairall 9. This value is almost one order of magnitude 
lower than the cooling rate derived for the one-tenth- 
solar metallicity HVC Complex C (log Ic ~ —27), but is 
relatively close to the values derived for solar-metallicity 
clouds in the Milky Way disk and in the lower Galactic 
halo (log Ic ~ —26; see Lehner, Wakker & Savage 2004, 
their Table 4). This result underlines the importance of 
metal cooling for the thermal state of metal-rich circum- 
galactic gas absorbers in the local Universe. 

5. DISCUSSION 

5.1. Enrichment history of the Magellanic Stream 

5.1.1. Alpha and nitrogen abundances 

One major result of our study is the surprisingly high 
metallicity of the MS toward Fairall 9. From the mea- 
sured sulfur abundance of [S/H] = -0.30 ± 0.04 fol- 
lows that the a abundance in the Stream in this di- 
rection (/ = 295, b = —58) is as high as 0.5 solar, 
which is ^ 5 times higher than the a abundance derived 
for other sightlines passing the MS toward NGC 7469 
{I ^ 347,6 = -64), RBS144 {I = 299, & = -66), and 
NGC 7714 (l = 88,6 = -56; Fox et al. 2010, 2013). Our 
detailed analysis of the ionization conditions in the gas, 
the comparison between Hi 21cm measurements from 
different radio telescopes with different beam sizes, and 
the analysis of the H i Ly a absorption toward Fairall 9 
(see Appendix) do not provide any evidence, that this 
high sulfur abundance could be a result from the various 
systematic errors that come along with our analysis (see 
also the discussion in Gibson et al. 2000 on this topic). 
Wc thus are forced to conclude that the measured high 
sulfur abundance in the gas reflects the true chemical 
composition of the Stream in this direction. 

As a guide to the following discussion, we have plot- 
ted in Fig. 9 the derived MS abundances (black filled 
circles) together with the SMC and LMC present-day 
stellar and nebular abundances (Russel & Dopita 1992; 
Hughes et al. 1998; red (SMC) and green (LMC) filled 
circles). With a sulfur abundance of 0.5 solar, the MS 
toward Fairall 9 exhibts an a abundance that is higher 
than the average present-day a abundance in both SMC 
and LMC (~ 0.3 and ~ 0.4 solar, respectively; Russell & 
Dopita 1992). This finding has profound implications for 
our understanding of the origin of the Magellanic Stream 
and its enrichment history. 

Early tidal models suggest that the MS was stripped 
from the Magellanic Clouds ^ 2 Gyr ago (Gardiner & 
Noguchi 1996). More recent models, that take into ac- 
count that the Magellanic Clouds possibly are on their 



first pass through the Milky Way halo, have included 
other relevant mechanisms that would help to unbind 
the gas, such as blowouts and tidal resonances (Con- 
nors et al.2006; Nidever et al.2008; Besla et al.2010). 
One would expect that the Stream's chemical composi- 
tion must reflect that of its parent galaxy at the time, 
when the gas was separated from the Magellanic Clouds, 
since the MS does not contain stars that would increase 
its metallicity since then. Following the well-defined age- 
metallicity relations of SMC and LMC (Pagel & Taut- 
vaisiene 1998; Harris & Zaritsky 2004; 2009), one would 
expect that the Stream's present-day a abundance is at 
most ^ 0.2 solar if originating in the SMC and ^ 0.25 
solar if originating in the LMC. While the observed MS 
abundances toward RBS 144, NGC 7469, and NGC 7714 
are in line with these abundance limits (and actually 
favour an SMC origin for this part of the Stream; see Pa- 
per I), the sulfur abundance in the MS toward Fairall 9 
is at least twice as high as expected from a simple tidal 
model that assumes a homogeneous pre-enrichment of 
the SMC/LMC gas (e.g., Dufour 1975) before the Stream 
was stripped off. 

A second important finding of our study, that adds to 
this puzzle, is the relatively low nitrogen abundance in 
the MS toward FairaU9. The [N/S] ratio (= [N/a] ratio) 
is —0.85 dex, which is very low for the relatively high 
a abundance of [q;/H]= —0.3 when compared to nearby 
extragalactic H ii regions and high-redshift damped Ly- 
man a (DLA) systems (see Pettini et al. 2008, Jenkins 
et al. 2005, and references therein). The a-process el- 
ements O, S, and Si are believed to be produced by 
Type II supernovae from massive progenitor stars, while 
the production of N as part of the CNO cycle in stars 
of different masses is less simple (and not yet fully un- 
derstood). The so-called "primary" nitrogen production 
occurs, when the seed elements C and O are produced 
within the star during the helium burning phase, while 
the for "secondary" nitrogen production these seed ele- 
ments were already present when the star condensed out 
of the ISM (Pettini ct al. 2002; Henry & Prochaska 2007). 
Primary N is believed to be produced predominantly 
by stars of intermediate masses on the asymtotic giant 
branch (e.g., Henry et al. 2000). Because of the longer 
lifetime of intermediate-mass stars, primary N therefore 
is expected to be released into the surrounding ISM with 
a time delay of ~ 250 Myr compared to a elements. At 
low metallicities less than ~ 0.5 solar (i.e., at LMC/SMC 
metallicitics ^1 — 2 Gyr ago, when the MS was sepa- 
rated from its parent galaxy) the nitrogen production is 
predominantly primary (see Pettini et al.2008). 

The low [N/S] and high [S/H] ratios observed in the MS 
toward Fairall 9 therefore indicate an abundance pattern 
that is dominated by the a enrichment from massive stars 
and Type II supernovae, while only very little (primary) 
nitrogen was deposited into the gas. 

5.1.2. Origin of the gas 

The most plausible scenario that would explain such 
an enrichment history, is, that the gas that later became 
the Stream was locally enriched in the Magellanic Clouds 
^1 — 2 Gyr ago with a elements by several supernova 
explosions in a star cluster or OB association, and then 
separated/stripped from the stellar body of the parent 
galaxy before the primary nitrogen was dumped into the 
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Fig. 9. — Comparison of present-day metal abundances in the 
SMC (blue dots), LMC (green dots), and in the Magellanic Stream 
toward Fairall9 (black dots). For SMC and LMC the metallicitics 
are derived from stellar and nebular abundances (from Russell & 
Dopita 1992; Hughes et al. 1998). 

gas and the metals could mix into the ambient interstel- 
lar gas. Eventually, the supernovae explosions may have 
pushed away the enriched material from the stellar disk 
so that the gas was already less gravitationally bound at 
the time its was stripped an incorporated into the Mag- 
ellanic Stream. In this scenario, it is the present-day N 
abundance in the MS that defines the metallicity floor of 
the parent galaxy and thus provides clues to the origin 
of the gas. 

Chemical evolution models of the Magellanic Clouds 
suggest that the mean metallicitics of SMC and LMC 
1 — 2 Gyr ago were 0.2 — 0.3 dex lower than today (Pagel 
& Tautvaisiene 1998; Harris & Zaritsky 2004, 2009). 
Therefore, and because of the substantially higher mean 
N abundance in the LMC compared to the SMC (see 
Fig. 9), the low N abundance measured in the Stream 
toward Fairall9 favours a SMC origin of the gas. An 
LMC origin cannot be ruled out, however, because the 
scatter in the present-day N abundances in the LMC 
is large (see Sect. 4.2) and the age-metallicity relations 
from Pagel & Tautvaisiene (1998) and Harris & Zarit- 
sky (2004,2009) may not apply to nitrogen because of its 
special production mechanism. In fact, considerung the 
position of the Fairall9 sightline (see Paper I, Fig. 1) and 
the radial velocity of the MS absorption in this direction, 
the bifurcation model of Nidever et al. (2008) predicts 
that the high- velocity gasous material towards Fairall 9 is 
part of the LMC filament of the Stream und thus should 
have a different chemical composition than the SMC fil- 
ament traced by the other MS sightlines (Paper I). An- 
other aspect that may be of relevance in this context is 
the fact that the Fairall 9 sightline lies only 14.3° on the 
sky from the SMC. The relative high metal-abundance in 
the Stream toward Fairall 9 thus may refiect the increas- 
ing importance of continuous ram-pressure stripping of 
metal-enriched SMC gas as the MCs get closer to the 
Milky Way. 

While the above outlined enrichment scenario appears 
feasible to explain qualitatively the observed trend seen 
in the abundance pattern in the MS toward Fairall 9, 
the question arises, how many massive stars and subse- 



quent supernova explosions would then be required to lift 
the [S/N] and [S/H] ratios to such a high level. To an- 
swer this question, one would first need to have an idea 
about volume and mass of the gas that is enriched in 
this manner. The fact that Fox et al. (2010, 2013) have 
determined a much lower a abundance in the Stream of 
'^0.1 solar toward three other MS sightlines using simi- 
lar data possibly suggests, that the high metallicity and 
low [N/S] ratio toward Fairall 9 represents a relatively 
local phenomenon. 

In Sect. 4.3, we have estimated a thickness of the ab- 
sorbing neutral gas layer in the Stream of d ~ 30 pc 
from the mean gas density (njj = 1-0 cm~^) and the 
total column density, assuming that the dust depletion 
values of Ti and Ca are similiar in the gas. In the fol- 
lowing, we regard this value as a realistic lower limit for 
the true physical size of the absorbing gas region. If we 
adopt the SN yields from Kobayashi et al. (2006) and as- 
sume a spherical symmetry for the absorbing gas region 
with a diameter and density as given above, we calculate 
that only a handful of massive stars would be required 
to boost the a abundance from initially 0.1 solar to 0.5 
solar in such a small volume of gas. A much higher num- 
ber of high-mass stars would be required, however, if the 
volume to be enriched would be substantially larger, as 
we outline in the following. Because the total mass of 
a spherical gas cloud can be written as the product of 
cloud volume and gas density (M oc Vnu oc d^nn), and 
because we use the relation d = Nn/nn to estimate d, 
we can write M oc n^ for a fixed (measured) hydrogen 
column density in the gas. Therefore, if the mean gas 
density in the MS toward Fairall 9 would be 0.1 cm~'^ 
rather than the assumed 1.0 cm~^, then the diameter of 
the cloud would be 300 pc instead of 30 pc, but it would 
require 100 times more massive stars to enrich this large 
volume to the level required. In view of the stellar con- 
tent of the most massive star forming regions in the Mag- 
ellanic Clouds (e.g., 30Doradus; Melnick 1985) and the 
possible presence of a major star burst in the MCs at the 
time, when the Magellanic Stream was formed (Harris & 
Zaritsky 2009; Weisz et al. 2013), a size of a few dozen up 
to a few hundred pc together with the proposed enrich- 
ment scenario is fully consistent with our understanding 
of the star-formation history of the Magellanic Clouds. 

In summary, the observed abundance pattern in the 
MS toward Fairall 9 suggests the presence of a high- 
metallicity gas filament in the Stream in this direc- 
tion, possibly originating in a region with enhanced star- 
formation activity in the SMC or LMC ~ 1 - 2 Gyr 
ago. Our study indicates that the chemical compostion 
and thus the enrichment history of the MS seems to be 
more complex than previously thought. The Stream pos- 
sibly is composed of different gas layers that have differ- 
ent chemical compositions and that originate from dif- 
ferent regions with local star-formation histories within 
the Magellanic Clouds. 

Interestingly, the MS absorber toward Fairall 9 is not 
the only example of a circumgalactic gas structure in the 
nearby Universe that exhibits a high overall metallicity 
together with a low N/a abundance ratio. For exam- 
ple, Jenkins et al. (2005) have measured [a/H]= —0.19 
and [N/a] < 0.59 in an intervening Lyman-limit system 
(LLS) at z = 0.081 toward the quasar PHL1811. This 
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absorption system appears to be associated with two 
neary spiral galaxies at impact parameters p = 34 h^Q 
and /) = 87 hj^ kpc. It possibly represents gaseous 
material that has been ejected or stripped from these 
galaxies (or other galaxies nearby) and thus possibly has 
an origin that is very similar to that of the Magellanic 
Stream. Note that low [N/a] ratios are also found in 
low-metallicity HVCs, such as Complex C (e.g., Richter 
at al. 2001b). Clearly, a systematic study of N/a ratios in 
circumgalactic metal absorbers at low z with HST/COS 
would be an important project to investigate, whether 
similar abundance patterns are typical for circumgalac- 
tic gaseous structures in the local Universe. 

5.2. Physical conditions and small-scale structure in the 

gas 

In addition to the very interesting chemical properties, 
the combined UVES/COS/Ff/S'.B/GASS data set pro- 
vides a deep insight into the physical conditions in the 
Stream toward Fairall 9. The data show a complex mul- 
tiphase gas structure that possibly spans a large range 
in gas densities and ionization conditions. 

The presence of H2 absorption toward Fairall 9 and 
NGC3783 (Sembach et al.2001) indicate that the MS 
and its Leading Arm hosts a widespread (because of the 
large absorption cross section/detection rate) cold neu- 
tral gas phase, possibly structured in a large amount of 
small, dense clumps or filaments. Assuming a total area 
of the MS and its Leading Arm of ~ 1500 deg^, and a 
typical size for these clumps of a few pc, there could be 
millions of these dense structures in the neutral gas body 
of the Stream, if the two H2 detections toward Fairall 9 
and NGC 3783 reflects the true absorption-cross section 
of this gas phase throughout the Stream's neutral gas 
body. 

The fact that the MS can maintain significant amounts 
of II2 at moderate gas densities (tt-h = 1 — 10 cm^'^) 
probably is a result of the relatively low intensity of 
the dissociating UV radiation at the Stream's location 
due to the absence of local UV sources (see also Fox et 
al. 2005). Note that the low excitation temperature of 
II2 for J > 2 (see also Sembach et al. 2001) provides in- 
dependent evidence for a relatively weak ambient UV ra- 
diation field, and thus supports this scenario. In the MS, 
the low dissociation rate is expected to overcompensate 
the low formation rate (because of the lower metallic- 
ity and dust content compared to the Milky Way ISM), 
thus resulting in molecular gas fractions that are rela- 
tively high for the given total gas column. With a grain- 
formation rate of Rms = 3 x 10~^^ cm^ s~^ and a density 
of uh = 5 cm~^ (see Sect. 4.4), the H2 formation time 
is long, fform — (^ms"h)~^ ~ 2 Gyr, which matches 
the tidal age of the Stream. Therefore, it seems most 
likely, that the II2 (or at least some fraction of it) has 
already formed in the parent galaxy and then survived 
the subsequent tidal stripping. A similar conclusion was 
drawn from the II2 observation in the Leading Arm of 
the MS toward NGC 3783 (Sembach et al.2001). Elab- 
orating our idea, that the gas was locally enriched by 
massive stars shortly before the Stream was separated 
from its parent galaxy, one could imagine the presence 
of relatively dense, compressed shells and fragments that 
formed the seed structures for the formation of II2 that 



then were carried along with the Stream's gaseous body. 
Independent constraints on the physical conditions in 
the gas come from the observed Nai/Caii ratio. We mea- 
sure a Nai/Caii column-density ratio of 0.23 in compo- 
nent 4, where the II2 is expected to reside. In the Milky 
Way ISM, such a ratio is typical for the warm neutral 
medium (WNM), where nn < 10 cm^^ a,nd T = 10^-10'' 
K, typcially, and Call serves a trace species (e.g., Craw- 
ford 1992). In such gas, and without dust depletion of 
Ca and Na, a nearly constant Nai/Caii of 0.025 ra- 
tio would be expected from detailed ionzation models of 
these ions (see Crawford 1992; Welty et al. 1996; Richter 
et al. 2011). If these numbers would apply also to the 
conditions in the MS, they would indicate relative dust 
depletions of Ca and Na of log 6{Ca.)— log 5(Na) « 0.9, 
or log (5(Ca) « 0.9 for log (5(Na) — 0. These values are in 
very good agreement with the dust-depletion estimates 
from our Cloudy model discussed in Sect. 4.3. 

5.3. Relevance to intervening QSO absorbers 

The Magellanic Stream represents a prime example for 
a high-column density circumgalactic tidal gas stream 
in the local Universe. If the Stream would be seen as 
QSO absorber, it would be classified as LLS, sub-DLA 
(sub-Damped Lyman a absorber), or DLA, depending 
on the position of the sightline passing through the gas. 
The H I column-density map presented in Fig. 1 provides 
an estimate on the covering fractions of these column- 
density regime in the MS in the direction of Fairall 9. The 
results from our multi-sightline campaign to study the 
properties of the MS therefore are also of relevance for the 
interpretation of interevening absorption-line systems at 
low redshift. 

5.3.1. Inhomogeneity of absorbers 

The first important result from our study that is of rel- 
evance for QSO absorption-line studies is, that the phys- 
ical conditions and the chemical composition appear to 
vary substantially within the Stream. 

If such inhomogeneities were typical for tidal gas fila- 
ments around galaxies, then the interpretation of absorp- 
tion spectra from circumgalactic gas structures around 
more distant galaxies (in terms of metallicity, molecu- 
lar content, physical conditions, gas mass, origin, etc.) 
would be affiicted with large systematic uncertainties. In 
fact, most studies that aim at exploring the connection 
between galaxies and their surrouding circumgalactic gas 
are limited to single sightlines that pass the galaxy en- 
vironment at a random impact parameter (e.g., Thoni 
et al.2011; Tumlinson et al. 2011; Ribaudo et al.2011). 
Physical and chemical properties derived from singlc- 
sightline studies, however, may not be representative at 
all for the conditions in the general gaseous environment 
(in the same way as conditions in the MS along the 
Fairall 9 sightline are not representative for the Stream 
as a whole). The inhomogeneous chemical composition 
of the Stream implies that the metals possibly are not 
well mixed in the gas. Former abundance variations in 
the parent galaxy due to local star-formation events thus 
may have been frozen into the Stream's spatial metal 
distribution. This would be not too surprising, however, 
since the main process that stir up and mix the interstel- 
lar gas within galaxies (i.e., supernova explosions, stellar 



16 



P. Richter et al. 



winds) cannot take place in the Stream simply due to 
the lack of stars. 

5.3.2. Molecular absorption. 

The two detections of H2 absorption in the MS with 
FUSE toward FairaU9 and NGC 3783 (Sembach et 
al. 2001) indicate that tidal gas streams around galax- 
ies may typically host a widespread, cold gas phase that 
has a substantial absorption cross section. This scenario 
is supported by the recent detection of H2 absorption 
in another, more distant circumgalactic tidal gas stream 
beyond the Local Group (Crighton et al. 2013). These 
findings also remind us that the presence of H2 in an 
intervening absorber with a complex velocity structure 
and a high neutral gas column density does not neces- 
sarily mean that one traces a gaseous disk of a galaxy. 

The relatively large absorption cross section of H2 
in neutral gas structures around galaxies, as found in 
the Milky Way halo (see Richter 2006 and references 
therein), can be explained by the circumstance, that the 
physical conditions in these star-less gas clouds favour 
the formation of diffuse molecular structures even at rel- 
atively moderate gas densities. In particular, the relative 
low intensity of the local UV radiation field, the efficient 
process of H2 line self-shielding, and the lack of the dis- 
tructive influence of massive stars probably represent im- 
portant aspects in this context. 

H2 observations in circumgalactic gas clouds not only 
provide important information on the physical and chem- 
ical conditions in these structures, they also are of high 
relevance for our understanding of the physics of molecu- 
lar hydrogen in diffuse gas under conditions that are very 
differnt from that in the local ISM (see also Sembach 
et al. 2001). The transition from neutral to molecular 
gas, in particular, represents one of the most important 
processes that govern the evolutionary state of galaxies 
at low and high redshift. Detailed measurements of H2 
fractions and dust-depletion patterns in (star-less) tidal 
streams at low redshift can help to constrain the critical 
formation rate of molecular hydrogen in low-metallicity 
environments and thus could be of great importance to 
better understand the distribution and cross section of 
molecular gas in and around high-redshift galaxies. 

6. SUMMARY 

In this second paper of our ongoing series to study 
the Magellanic Stream in absorption we have analyzed 
newly obtained optical and UV absorption-line data 
from HST/COS and VLT/UVES together with archival 
FUSE and Hi 21cm emission data from GASS to study 
the chemical composition and the physical conditions 
in the Magellanic Stream in the direction of the quasar 
Fairall9. Our main results are the following: 

1. Metal absorption in the Magellanic Stream (MS) 
is detected in seven individual absorption components 
centered at wlsr = +143, -Kl63, -^182, -^194, +204 and 
+218 kms~^, indicating a complex internal velocity 
strucutre of the MS in this direction. Detected ions, 
atoms and molecules in the Stream include CiV, Siiv, 
Sim, Cii, Cii*, Alii, Sin, Sii, Can, Tin, Fen, Tin, 
Oi, Ni, Nai, and II2. 

2. From the unsaturated Sn absorption and a Cloudy 



photoionization model we measure an a abundance in 
the Stream of [S/H] = -0.30 ± 0.04 (0.50 solar), which 
is substantially higher than that found in the Stream 
along the lines of sight toward NGC 7469, RBS 144, 
and NGC 7714 (Fox et al.2010, 2013). Unfortunately, 
the unresolved, saturated Oi A1302.2 line cannot be 
used to independently constrain the a abundance 
in the MS toward Fairall9. Contrary to sulfur, we 
measure a very low nitrogen abundance in the gas of 
[N/H] = -1.15 + 0.06. The resulting [N/S] ratio is -0.85 
dex, which is very low compared to other environments 
with similarly high a abundances. The low [N/S] and 
high [S/H] ratios observed in the MS toward Fairall9 
suggest an abundance pattern that is dominated by 
the a enrichment from massive stars and Type II 
supernovae, while only very little primary nitrogen was 
deposited into the gas when the Stream was separated 
from the Magellanic Clouds. 

3. The detection of very narrow Nai and H2 absorption 
(with 6 « 2 kms^^) in the component at wlsr ~ +188 
kms~^ indicates the presence of a compact (pc-scale), 
cold gas stucure in the MS along this sightline. From 
the analysis of the newly reduced archival FUSE data 
of Fairall9 we measure a total molecular hydrogen 
column density of log 7V(H2) ~ 17.93lg ijg, improving 
previous results from Richter et al. (2001a). From the 
analysis of the H2 rotational excitation we obtain a 
kinetic temperature in the cold neutral gas phase of 
r « 93 K. For the gas density we derive n « 4 — 8 cm~^, 
assuming that the H2 gas is in a formation-dissociation 
equilibrium. The resulting estimate for the thermal gas 
pressure is P/k ~ 750 cm^'^ K, thus in good agreement 
with values derived from previous studies of the Stream. 
The detection of H2 absorbing structures in the MS, 
whose linear and angular sizes must be very small 
(~ 1 pc and ^ 3 arcseconds), indicates that the neutral 
gas body of the Stream is highly clumped and structured. 

4. The relatively low column densities of Fen, Ti,n, 
Nin, Aln, Ca,n indicate that the gas phase abundances 
of these elements are affected by dust depletion. We 
combine our column-density measurements for these ions 
with a Cloudy photoionization model and derive dust 
depletion values relative to sulfur of log (5(Si) — —0.27, 
log (5(Fe) == -0.56, log 5{Ki) < -0.51, log S{T[) = -1.32, 
log (5(A1) = -0.62, and log (5(Ca) = to -1.58. These 
depiction values are similar to those found in warm, 
diffuse clouds in the lower Milky Way halo (Savage & 
Sembach 1996). 

5. Our study indicates that the enrichment history of 
the Magellanic Stream as well as the physical conditions 
in the Stream are more complex than previously known. 
The abundances and gas-to-dust ratios measured in 
the Strcan along the Fairall9 sightline are substantially 
higher than what is found along other MS sightlines. 
The low N abundance favours an SMC origin of the MS 
gas toward Fairall9, although an LMC origin cannot 
be ruled out. The high sulfur abundance in the gas 
possibly indicates a massive a enrichment from massive 
stars in a region of enhanced star-formation ^1 — 2 Gyr 
ago and then was stripped from the Magellanic Clouds 
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and incorporated into the MS, before the gas could be 
enriched in nitrogen from intermediate stars. 
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Fig. 10. — Hi Lyo absorption in the COS spectrum of Fairall9 in the wavelength range between 1195 and 1235 A. The red solid line 
indicates a two-component absorption model with a Milky Way Hi column density of log Af(Hl) = 20.25 and a MS Hi column density of 
log Ar(Hi) = 19.95. 



APPENDIX 
HI LYa ABSORPTION 

In Fig. 10 we show the normalized COS spectrum in the wavelength range between 1195 and 1235 A. Strong Hi 
Lya absorption centered at 1216 A from neutral gas in the Milky Way (MW) and in the Magellanic Stream with 
well-defined damping wings is visible. The 21cm data from various observations and instruments (Sect. 3.3; Gibson et 
al. 2000) indicate that the Milky Way disk gas dominates the neutral-gas column density along this sightline 

The strength and shape of the Ly a absorption toward Fairall 9 supports this conclusion, but the relative contribution 
of MW and MS cannot be tightly constrained from a Voigt-profile fit of the Ly a absorption. The red solid line in 
Fig. 10 shows a two-component absorption model with a MW Hi column density of log iV(Hi) = 20.25 and a MS Hi 
column density of log iV(Hi) — 19.95. Therefore, the Lya absorption profile is fully consistent with the value of log 
iV(Hi) = 19.95 in MS, as indicated by the GASS 21cm data. However, the Lya profile is also consistent with MS 
column densities in the range log 7V(Hi) = 19.75 — 20.15, if the MW contribution is modified accordingly. In conclusion, 
the Lya absorption does not provide constraining limits on the Hi column density in the Magellanic Stream toward 
Fairall 9 that could help to minimize the systematic error for the determination of log A^(Hi), as discussed in Sect. 3.3. 

ABSORPTION MODELING OF UV LINES 

In Sect. 2.5 we have outlined the method that we have used to model the UV fines of the low ions in the COS spectrum 
of Fairall 9 to infer column densities of these species. We here present additional information on the modeling procedure 
to provide an insight into the allowed parameter range for (Ni, hi) and A^tot for each ion. 

In Fig. 11 we show velocity plots of the low ions from synthetic spectra that have been generated by us based on 
different input models, in which we have systematically varied the parameters Ni and hi. The resulting model profiles 
are compared to the best-fitting model presented in Fig. 3 and Table 1. 

In the left panel of Fig. 11 (a), the best-fiting model is indicated with the black solid line, while the lines in red, 
green, and blue show models, in which all column densities Ni in the individual subcomponents have been increased 
by 0.3, 0.5, and 1.0 dex, respectively. The velocity structure and the b values remain unchanged. In the middle panel 
of Fig. 11 (b), the column densities Ni in the individual subcomponents have been decreased by 0.3, 0.5, and 1.0 dex 
(red, green, blue), again without changes in Vi and hi. These two figures indicate that the shapes of the strong and 
saturated lines from Cii, Oi and Sin are basically not affected by the change in column density. In contrast, the 
absorption depths in the weaker lines of Ni, Alii, Sill, Sii, and Fen are strongly affected by the increase or decrease 
of the column densities with respect to the reference model (black) , demonstrating that these lines represent sensitive 
diagnostics to constrain A^^ and A'tot for these ions in the Stream. 

In the right panel of Fig. 11 (c) we show a model in which we have increased the Doppler parameters hi for the 
individual subcomponents by factors of 1.2 (red), 1.5 (green), and 2.0 (blue) with respect to the reference model 
(black), while the column densities and the overall velocity structure remain unchanged. Because of the relatively low 
spectral resolution of COS compared to the absolute values of hi, the changes in ht are irrelevant for the shape of the 
weak, unsaturated lines, while the strong, saturated lines mildly grow in strength in their absorption wings. Therefore, 
even if we would underestimate the values of hi for the singly-ionized species Aln, Sin, Sil, and Fell when adopting 
the values from Can, as one could argue because these ions may trace a slightly different (more extended) gas phase, 
our column-density estimate would not be affected at all by this systematic error. 

These figures, together with Fig. 2, underline, that the total gas column densities of Ni, Aln, Sin, Sii, and Fen in 
the MS, as listed in Table 1, are well constrained by the COS data and our absorption model. 
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Fig. 11. — Velocity profiles of low ions, as derived from synthetic spectra from different input models for the gas absorption in the 
Magellanic Stream toward Fairall 9. Details are given in the text in the Appendix. 



